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The Dyson— Schwinger equations of QCD, truncated to ladder-rainbow level, are 
used to calculate meson form factors in impulse approximation. The infrared 
strength of the ladder-rainbow kernel is described by two parameters fitted to the 
chiral condensate and /•„■; the ultraviolet behavior is fixed by the QCD running 
coupling. This obtained elastic form factors F 7r (Q 2 ) and Fk(Q 2 ) agree well with 
the available data. We also calculate the p — > iry and K* — > transition form 
factors, which are useful for meson-exchange models. 



1. Dyson Schwinger equations 

The set of Dyson-Schwinger equations [DSEs] form a useful tool to obtain 
a microscopic description of hadronic properties 1 . Here we use the DSEs 
to calculate elastic and transition form factors of the light mesons. The 
dressed quark propagator, as obtained from its DSE, together with the 
meson Bethe-Salpeter amplitude [BSA] and the dressed quark-photon ver- 
tex, form the necessary elements for calculations of form factors in impulse 
approximation, such as the pion elastic form factor 2 . 

The DSE for the renormalized quark propagator in Euclidean space is 1 

-JLtfD^lp-ql^pSb) Tl(q,p) ,(1) 

where D f _ ll/ (k) is the dressed-gluon propagator and T\,(q;p) the dressed- 
quark-gluon vertex. The most general solution of Eq. (1) has the form 
S(p)^ 1 — i^>A(p 2 ) + B(p 2 ) and is renormalized at spacelike /i 2 according to 
A(/i 2 ) = 1 and B(p 2 ) = m(/x) with m(//) the current quark mass. 

Mesons are described by solutions of the Bethe-Salpeter equation [BSE] 

■j-?- K(p, q; Q) S(q+) T H (q + , g_; Q) S{q.) , (2) 

at discrete values of Q 2 = —m 2 Hl where mjy is the meson mass. In this 
equation, p + = p + rjQ and p_ = p — (1 — rj)Q are the outgoing and 
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incoming quark momenta respectively, and similarly for q±. The kernel K 
is the renormalized, amputated qq scattering kernel that is irreducible with 
respect to a pair of qq lines. Together with the canonical normalization 
condition for qq bound states, Eq. (2) completely determines the bound 
state BSA Different types of mesons, such as pseudoscalar or vector 
mesons, are characterized by different Dirac structures. 

The quark- photon vertex, r M (p+,p_; Q), with Q the photon momentum 
and p± the quark momenta, is the solution of the inhomogeneous BSE 



Solutions of the homogeneous version of Eq. (3) define vector meson 
bound states at timelike photon momenta Q 2 = —my. It follows that 
T^(p + , p_; Q) has poles at these locations 2 ' 3 . 
To solve the BSE, we use a ladder truncation, 



in conjunction with the rainbow truncation for the quark DSE, Eq. (1): 



K(q,P) -> 7.AV2 and Z l9 2 D^(k) - 4na cS (k 2 )D%(k). Here, D%{k) is 



the free gluon propagator in Landau gauge, and a eS (k 2 ) the effective quark- 
quark interaction, which reduces to the one-loop QCD running coupling 
a 1_loop (fc 2 ) in the perturbative region. For the infrared behavior of the 
interaction, we employ an Ansatz 4 ' 5 that is sufficiently strong to produce 
a realistic value for the chiral condensate of about (240 GeV) 3 . The model 
parameters 5 , along with the quark masses, are fitted to give a good descrip- 
tion of the chiral condensate, m^/x and f n . The obtained results for the 
light vector meson masses are within 5% of their experimental values, and 
the vector meson electroweak decay constants are within 9% of the data 5 . 

This truncation preserves both the vector Ward-Takahashi identity 
[WTI] for the qq^y vertex and the axial-vector WTI, independent of the 
details of the effective interaction. The latter ensures the existence of 
massless pseudoscalar mesons associated with dynamical chiral symmetry 
breaking 4,6 . In combination with impulse approximation, the former en- 
sures electromagnetic current conservation 3 . 

2. Meson Form Factors 

In impulse approximation, meson form factors are generically described by 



F»(p+,P-;Q) = ^2 7m + 




K(p, q; P) - -a cS ((p - q) 2 ) D%{p - q)^ ® f 7i 



(4) 



I abc (P,Q,K) =N C 




Tr[S a (q)T ab (q, q ';P)S b (q') 



(27T) 4 
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where q — q' = P, q' — q" = Q, q" — q = K , and momentum conservation 
dictates P + Q + K = 0. In Eq. (5), S l is the dressed quark propagator 
with flavor index i, and r y (fc, k';P) stands for a generic vertex function 
with incoming quark flavor j and momentum fc', and outgoing quark flavor 
i and momentum k. Depending on the specific process under consideration, 
this vertex function could be a meson BSA or a quark-photon vertex. In 
the calculations discussed below, the propagators and the vertices are all 
obtained as solutions of their respective DSE in rainbow-ladder truncation, 
without adjusting any of the model parameters. 



2.1. Pion and kaon elastic form factors 

There are two diagrams that contribute to meson electromagnetic form 
factors: one with the photon coupled to the quark and one with the photon 
coupled to the antiquark respectively. With photon momentum Q, and 
incoming and outgoing meson momenta P =F Q/2, wc can define a form 
factor for each of these diagrams 3 

2P u F all {Q 2 ) =lf b (P-Q/2,Q,-(P + Qm)- (6) 

We work in the isospin symmetry limit, and thus F 7T (Q 2 ) = F U u U {Q 2 )- 
The K + and K° form factors are given by F K + = ^F uSu + and 
F K o = -\F dsd + ±F dsg respectively. 

Our result for Q 2 F 7T and F K + are shown in Fig. 1, together with the 
experimental data 7 ' 8 ' 9 ' 10 . Up to about Q 2 = 2GeV 2 , our result for F^ can 
be described very well by a monopole with our calculated p-mass, m p — 
742 MeV (note that our calculated p-mass is slightly below the experimental 




1 



value). Above this value, our curve starts to deviate more and more from 
this naive VMD monopole. Our result is in excellent agreement with the 
most recent JLab data 9 ; it would be very interesting to compare with future 
JLab data in the 3 to 5 GeV 2 range, where we expect to see a significant 
deviation from the naive monopole behavior. 

Also our results for Fr agree with the available experimental data 10 , as 
do both the neutral and the charged kaon charge radius 3 . Again, our curve 
for Fk can be fitted quite well by a monopole up to about Q 2 = 2 GeV 2 , see 
Fig. 1; above 2 GeV 2 , we see a clear deviation from a monopole behavior. 
Also for this form factor it would be interesting to compare with future 
JLab data at larger Q 2 . 

2.2. Vector-pseudoscalar-photon transitions 

We can describe the radiative decay of the vector mesons using the same 
loop integral, Eq. (5), this time with one vector meson BSA, one pseu- 
doscalar BSA, and one 997- vertex 11 . The on-shell value gives us the cou- 
pling constant. For virtual photons, we can define a transition form factor 
Fyp-yiQ 2 ), normalized to 1 at Q 2 = 0, which can be used in estimating 
meson-exchange contributions to hadronic processes 12 ' 13 . 

In the isospin limit, both the p° tt° 7 and p ± 7r ± 7 vertices are given by 

\l%r(P,Q,-(P + Q)) = ^ e^ aP P a QpF p ^{Q 2 ), (7) 

where P is the p momentum. The ujitj vertex is a factor of 3 larger, due 
to the difference in isospin factors. For the K* — ► decay, we have to 
add two terms: one with the photon coupled to the s-quark and one with 
the photon coupled to the u- or <i-quark, corresponding to the charged or 
neutral K* decay respectively 11 . 

As Eq. (7) shows, it is gvPj/ m v that is the natural outcome of our cal- 
culations; therefore, it is this combination that we give in Table 1, together 
with the corresponding partial decay widths 11 . The agreement between 
theory and experiment for gyp^jmv is within about 10%, except for the 
discrepancy in the charged K* — ► if 7 decay for which we have no explana- 
tion. Likewise the large difference between the neutral and charged p decay 
width is beyond the reach of the isospin symmetric impulse approximation. 



Table 1. Vector meson radiative decays: g/m in GeV 1 and Yy^p 1 in keV. 





g/m 




g/m 




g/m 




g/m 


r K*°K<>-/ 


calc. 


0.69 


53 


2.07 


479 


0.99 


90 


1.19 


130 


expt. 14 


0.74 


68 


2.31 


717 


0.83 


50.3 


1.28 


116 
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Note that part of the difference between the experimental and calculated 
decay width conies from the phase space factor because our calculated vec- 
tor meson masses deviate up to 5% from the physical masses. 

The corresponding transition form factors are shown in Fig. 2. In con- 
trast to the elastic form factors, these transition form factors fall off sig- 
nificantly faster than a VMD-like monopole; our numerical results 11 for 
FpTT-yiQ 2 ) suggest an asymptotic behavior of 1/Q 4 . Only in the timelike 
region, near the vector meson pole, do we see a clear VMD-like behavior. 

For the K*K^f form factor the situation is more complicated due to 
the interference of the diagrams with the photon coupled to the s-quark 
and to the u- or d-quark. This causes the charged form factor to fall off 
much more rapidly than both the neutral K*K^y form factor and the elastic 
form factor Fk(Q 2 ), as can be seen from Fig. 2. The latter implies that 
the contribution to charged kaon electroproduction from intermediate K* 
exchange gets suppressed with increasing Q 2 compared to the contribution 
from virtual kaon exchange. 



2.3. Remarks on meson electroproduction 

We can use this approach to estimate the range of validity of meson- 
dominance models 11 . For off-shell pions, the meson-dominance assumption 
appears to be quite good for spacelike momentum transfer of the order 
of t ~ 0.1 GeV 2 . On the other hand, for heavier mesons such as kaons 
or p-mesons, the naive meson-dominance assumption introduces significant 
errors, even at small spacelike values of t <~ 0.2 GeV 2 . In addition, any 
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meson-exchange model to describe meson electroproduction necessarily in- 
troduces off-shell ambiguities. Clearly, a microscopic description is needed. 

Any microscopic description of meson-electroproduction requires a 
quark-gluon description of the nuclcon. It is a difficult task to combine 
such a description with the meson form factors considered here. Recently, 
significant progress in describing processes involving four external particles 
has been made, using the rainbow-ladder truncation of the DSEs, in con- 
junction with an extention of the impulse approximation 16 . This approach 
incorporates the non-analytic effects of intermediate meson- exchange con- 
tributions while avoiding ambiguous off-shell definitions. It has succesfully 
been applied to ir-ir scattering, where we can identify the non-analytic 
contributions from a- and p-exchange to the S-wave and P-wave scattering 
amplitudes respectively 17 . Furthermore, it was shown to reproduce the cor- 
rect chiral limit 16 . We plan to extend this approach to describe processes 
such as 77r7T7r, pion compton scattering, and meson electroproduction. 
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